Extreme fluctuations and the flnite lifetime of the turbulent state 
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We argue that the transition to turbulence is controlled by large amplitude events that follow 
extreme distribution theory. The theory suggests an explanation for recent observations of the 
turbulent state lifetime which exhibit super-exponential scaling behaviour with Reynolds number. 
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The fundamental nature and stability of the turbulent 
state of fluids remains an open and challenging question. 
Fluid flow is characterized by a dimensionless number Re, 
which depends on the characteristic length L, velocity U 
and kinematic viscosity of the fluid v through the rela- 
tion Re = UL/v. As the Reynolds number increases from 
zero, the flow becomes increasingly structured and even- 
tually statistical in nature, and at large Re, the flow is 
said to be turbulent The conventional assumption — 
that the turbulent state is absolutely stable—has been 
challenged recently by a series of theoretical 12] and ex- 
perimental probesjSHTl of the transition to turbulence. 
Taken as a whole, these works suggest that turbulence 
might, in some flow reg imes at least, be a long-lived 
metastable stated, Such a view would be con- 

sistent with the fact that long-lived transient turbulent 
states can be excited as finite-amplitude instabilities of 
the laminar state, so that the laminar and turbulent 
states can coexist (for a review of foundational work in 
this area, see e.g. Ref. Ildh recent developments are 
summarized in Refs. [9l-lllL Il3l|). However, the question 
remains as to whether the turbulent state is ever sustain- 
able with an infinite lifetime for finite Reynolds numbers. 
This is a difficult experimental question to decide, be- 
cause the lifetime of the turbulent state can become so 
long that measurements become impossible. With the 
necessary restriction to a small range of Reynolds num- 
bers, the data have, until recently, been difficult to inter- 
pret in a compelling way. 

In a set of elegant and remarkably accurate experi- 
ments on transitional pipe turbulence 6], Hof et al. have 
brought into question the idea that pipe fiow turbu- 
lence is stable at long times beyond a finite critical 
Reynolds number (l^.[l5|. The laminar state of a straight 
smooth pipe flow is linearly stable at all Reynolds num- 
bers (see e.g. Ref. (^16'])), but a sufficiently large per- 
turbation triggers localized turbulent puffs that persist 
for long times. The decay of the transient turbulent 
state is reported to follow a Poisson distribution, with 
a lifetime T(Re) that increases sharply with increasing 
Reynolds number. The new measurements of the life- 
time of these localized puffsQ reveal that T(Re) appar- 



ently only diverges at infinite Reynolds number, scaling 
in a super-exponential way with Re. Similar observations 
in another linearly stable fiow — Taylor-Couette flow with 
outer cylinder rotation — have recently been reported by 
Borrero-Echeverry et al. . 

In this article, we show that the form of the experi- 
mental data is consistent with a simple and general in- 
terpretation predicated on the use of extremal statis- 
tics. Our approach is related to the notion that the 
transient turbulent phenomena reflect escape from a low- 
dimensional dynamical attractor 17- 19], but we conceive 
turbulence as a spatially-extended phenomenon with a 
large number of degrees of freedom. The determining 
factor for the suppression of a puff is the probability that 
the largest fluctuation in a spatio-temporal interval con- 
sisting of multiple fluctuations fails to attain a threshold 
value. Thus we need to calculate the probability that 
the maximum amplitude of turbulent velocity fluctua- 
tions 6v{x,t) falls below some threshold value, which we 
term Be- We will assume below that once the turbulence 
has been sufficiently suppressed, the turbulent state is 
quenched, an assurngtion consistent with previously pub- 
lished analyses !20-22j. Our calculation shows that the 
super-exponential dependence of the lifetime of the tur- 
bulent state is a generic result of extremal statistics. 

In order to understand the lifetime of turbulent puffs, 
we assume that turbulent velocity configurations may be 
regarded as independent beyond a correlation time tq, 
and that there is a probability p that the puff will be 
suppressed within each time interval tq. Then, the life- 
time statistics will be Poisson. The probability P that 
turbulence persists to a time t after becoming established 
at a time Iq is P = (1 — p)*^, where the number of inter- 
vals is M — {t — to)/To. Therefore 

ln(P) -Mln(l-p) - l{t-to)\n{l-p), (1) 



To 



and so it follows that tq /t 
we can estimate ln(l — p) = 
lifetime in the form 



- In(l-p). Since 1 > p > 0, 
-p and therefore express the 



T = To/p 



(2) 
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where p depends on Re. 

We now determine how p depends on the Reynolds 
number of the flow and potentiaUy other factors. Within 
a spatial and temporal interval, multiple fluctuations oc- 
cur, sampled from the turbulent velocity distribution 
Pt{5v). The energy associated with these fluctuations 
is proportional to Sv^. We assume that when the energy 
fails to attain a certain threshold at all points in 

the puff, the turbulent state becomes unstable and de- 
cays. Thus if the largest velocity fluctuation is less than 
the threshold, all of the turbulent fluctuations are less 
than the threshold. Accordingly, it is necessary to calcu- 
late the probability distribution of the maximum of the 
fluctuations within a puff, in order to ensure that the 
largest fluctuation is below the threshold. 

We consider primarily a Gaussian distribution of ve- 
locity fluctuations, but our arguments also apply to the 
case of an exponential distribution (as is the case at high 
Reynolds numbers). We seek the probability distribution 
Pm{x) for the maximum x of a set of energy fluctuations 
{<5wf}, where i = \ . . . N . N represents the number of 
degrees of freedom, and should scale with the size of the 
turbulent puff, denoted here by A. Standard results from 
extreme statistics theory show that the appropriate result 
is the family of Fisher-Tippett distributions 23] . In par- 
ticular, the universality class for Pm must be the Type 
I Fisher-Tippett distribution, sometimes known as the 
Gumbel distribution [H, Hi] 

Pm{x) = ■^exp(-(x-/^)/^)exp(exp(-(x-Ai)//3)), (3) 

where /3 sets the scale and /j, the location of the distribu- 
tion. Note that the scale and location will depend on N , 
because the maximum of a set of random variables will 
be an increasing function of the number of random vari- 
ables. In particular, for the Gaussian case, Fisher and 
Tippett showed that asymptotically 



/i ~ VlogiV, (i ^ 1/VlogiV. 



(4) 



The mean and standard deviation of the Gumbel dis- 
tribution are /j, -I- r/3 and Ptt/^JQ respectively, where 
r sa 0.577 is the Euler-Mascheroni constant. The corre- 
sponding cumulative distribution is the probability that 
X < X, and is given by 

F{X) = f Pm{x) dx = exp(- exp(-(X-/i)//3)). (5) 

Thus, p — F{Bc), where Be is the threshold. 

We anticipate that Be is a decreasing function of Re, 
reflecting the intuition that at higher Re, turbulence can 
be more easily sustained by small fluctuations. We will 
consider the behavior of Be, as this sets the thresh- 
old in the distribution of energy maxima. The experi- 
ments are conducted in nominally smooth pipes within 



a narrow range of Re, so it is appropriate to expand 
Be about a particular Reynolds number Rcq, leading to 
Be = S^ + Bi(Re-Reo)-|-0(Re^), where B° and B^ are 
coeflicients. In order to describe the same Reynolds num- 
ber regime of the experiments, Rco may be interpreted to 
be a characteristic Reynolds number at which localized 
turbulent puffs first are observable, so that the lifetime 
is order tq. This onset is not a precisely defined point, 
but for concreteness we will specifically define it to be the 
Reynolds number where t = ctq. We will see that this 
choice simplifies the analysis below, but we emphasize 
that irrespective of whether or not we use the coefficient 
e or some other number of 0(1), our main predictions 
for the super-exponential distribution are not affected. 
The freedom of choice in this definition of Rcq is com- 
pletely analogous to the arbitrariness in the definition 
of the coexistence point between liquid and gas, which 
is also dominated by nucleation phenomena, as has also 
been noticed by Mannevillc[22]. 

Collecting results, we find that the average lifetime of a 
turbulent puff will have the approximate functional form: 

T = Toexp(exp(-(B^ + S^i(Re-Reo)-f 0((Re-Reo)^))) 

(6) 

in agreement with experimental findings. The coefficient 
To may in principle depend on pipe length or aspect ra- 
tio if these factors change the spatial scale on which re- 
gions of the localized turbulent puff are statistically in- 
dependent, and the timescale on which the state of the 
puff loses memory of previous states. The Taylor expan- 
sion only needs to be carried out to first order given the 
(understandably) small range of Reynolds number over 
which the experiments measuring the lifetime of turbu- 
lent puffs have been conducted. From ([5]) we can write 



lnln(r/To) = -{Be - fJ-)/P = ciRe -I- C2 



(7) 



where we have used the notation of Ref. ( 6]) to denote 
the coefficients ci and C2 of the linear fit to the data. 
Gomparing with Eq. ([6]) we read off that 

ci = ^Bl//3 C2 = -(SO-M-BiReo)//3 (8) 

Now, at Reo, the lifetime becomes comparable to the 
correlation time tq, and to be concrete, we chose that 
for Re = Rco, t = ctq, although it is straightforward to 
verify that our results have only a very weak dependence 
on the precise coefficient used. Then from Eq. ([6|), we 
see that Be = /i, and thus the ratio of the coefficients 
C2/C1 = — Rcq. The physical interpretation, if any, of 
the fltting parameter Reo is not clear to us, because al- 
though it is defined loosely as a characteristic Reynolds 
number below which the lifetime of the turbulent state is 
too small to be observable, any systematic dependence on 
intrinsic features of turbulence, flow geometry and per- 
haps wall roughness is beyond the scope of our work and 
of the experiments available at this time. Moreover, the 
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choice of definition of the Reynolds number in any partic- 
ular geometry is not unique once multiple length-scales 
are present, as in the case of Taylor-Couette flow, and 
thus it is hard to identify a unique and consistent defi- 
nition of Re and Rcq in order to retain meaning across 
different flow geometries. 

We conclude with a comment about the dependence of 
T on puff length A. The number of degrees of freedom 
N active in the turbulent puff is proportional to A. In 
our approach, the A-dependence of r can be estimated 
by substituting the scaling of /i and /3 as given by Eq. 
(|4]) into our formula for t. This yields that log(r) cx A to 
leading order, showing that the super-exponential scaling 
with Reynolds number does not translate into a super- 
exponential scaling with length, and is consistent with 
numerical measurements reported in Ref. ([21]). This 
prediction also applies if the probability distribution for 
velocity fluctuations is exponential, because in this case 
/i ^ logiV, i.e. as log A, but /3 does not scale with jv[23t. 

The interpretation of the lifetime statistics given here 
is related to that suggested recently by Manneville 221 , 
and earlier considerations by Pomeau on nucleation phe- 
nomena in hydrodynamics |2^. Indeed, any spatially- 
extended dynamical system with a memoryless sub- 
critical bifurcation should be expected to yield extremal 
statistics, and we hope to report on detailed calculations 
of this phenomenon in a future publication 27 [ . 
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